NOLTR  66-211 


CO 

X 

X* 

lO 

CD 

Q 

<! 


REPRESENTATIOIJ  OF  A  NONSPHERICAL 
UNDERWATER  EXPLOSION  BUBBLE  BY  A 
VORTEX  SHEET  AND  SOURCE  SHEET 


D  D  C 

ffj\P^r=nn  n-irj; 

JAN  2  6  1967  ‘ 

UucjLaU  criaily 

B 


n 


NOL 


22  NOVEMBER  1966 


UNITED  STATES  NAVAL  ORDNANCE  LABORATORY,  WHITE  OAK,  MARYLAND 


(NJ 

i 

vO 

X> 

Qd 


Distribution  of  this  document  is  unlimited. 


NOLTR  66-211 


REffiESElfEATI ON  OF  A  NONSPHERICAL  UNDERWATER  EXPLOSION 
BUBBLE  BY  A  VORTEX  SHEET  AND  SOURCE  SHEET 

by 

Larry  W.  Bell 

ABSTRACT:  The  equations  for  the  velocity  potential,  partial  time 
derivative  of  the  velocity  potential,  tangential  velocity,  and  normal 
velocity  for  a  source  sheet  and  vortex  sheet  representation  of  the 
nonspherical  underwater  explosion  bubble  have  been  derived  and  programmed 
for  the  IBM  79^0  computer.  These  equations  were  tested  for  the  known 
case  of  a  spherical  underwater  explosion  bubble  ana  the  resurts  were 
sufficiently  accurate  for  practical  application.  The  conclusion  is 
that  the  nonspherical  underwater  explosion  bubble  can  be  represented 
by  a  source  sheet  and  vortex  sheet.  This  representation  will  not  inter¬ 
fere  with  the  Jet  rising  into  the  interior  of  the  bubble,  as  is  the  case 
of  the  earlier  models  consisting  of  a  single  point  source  and  a  single 
dipole  located  on  the  axis  of  symmetry. 
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REPRESENTATION  OF  A  NON5PHERICAL  UNDERWATER  EXPLOSION 
BUBBLE  BY  A  VORTEX  SHEET  AND  SOURCE  SHEET 

INTRODUCTION 

After  the  detonation  of  an  explosive  charge  under  vater  and  the  emis¬ 
sion  of  the  shock  wave,  the  gaseous  products  of  the  explosion,  vhich  have 
f  jrmed  a  gas  bubble,  continue  to  expand  outward  at  a  gradually  decreasing 
rate.  The  result  of  this  expansion  is  a  radial  displacement  of  vater. 

When  the  radial  velocity  becomes  zero,  the  bubble  has  reached  its  first 
maximum  radius  and  the  pressure  in  the  bubble  is  far  below  the  hydrostatic 
pressure.  Subsequently,  the  bubble  contracts  and  the  bubble  pressure 
increases.  When  the  radial  velocity  again  becomes  zero,  the  bubble  pres¬ 
sure  is  substantially  greater  than  the  hydrostatic  pressure.  The  bubble 
has  now  reached  its  first  minimum  radius  and  begins  another  expansion 
phase.  The  gas  bubble  may  pulsate  in  this  manner  several  times  depending 
on  the  depth  at  Wiich  the  charge  is  detonated. 

Throughout  the  first  expansion  phase  of  the  bubble  there  is  very 
little  vertical  motion,  commonly  called  migration.  However,  near  the 
minimum  radius  an  appreciable  upward  displacement  may  occur.  The  migra¬ 
tion  is  a  result  of  the  buoyancy  of  the  bubble  and  depends  on  the  depth 
of  explosion  and  the  charge  weight.  There  is  little  migration  if  small 
charges  are  exploded  at  great  depths,  say  one  pound  at  200  feet.  However, 
the  explosion  of  large  charges  at  ordinary  depths,  charges  Of  the  order 
of  300  pounds  detonated  to  60  feet,  would  show  a  strong  migration. 

The  classical  treatment  of  the  underwater  explosion  bubble  is  based 
on  the  assumption  that  the  bubble  remains  spherical  at  all  times.  The 
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explosion  of  snail  charges  at  great  depths  essentially  behave  in  this 
manner.  However,  it  has  been  observed  in  the  explosion  of  large  charges 
at  ordinary  depths,  where  strong  migration  is  present,  that  the  bubble  does 
not  remain  spherical  but  undergoes  a  very  definite  change  in  shape. 
Initially  the  bubble  is  spherical  and  remains  spherical  sore what  beyond 
the  first  maximum  radius.  However,  as  the  bubble  contracts  to  its  first 
minimum,  the  bubble  is  flattened;  and  shortly  before  the  minimum,  the 
bottom  begins  to  move  upward  into  the  interior  of  the  bubble,  figure  1, 
and  eventually  impacts  with  the  upper  surface  cf  the  bubble,  figure  2. 

The  tongue  of  water  protruding  into  the  interior  of  th®  bubble  is  called 
the  bubble  Jet. 

Several  authors  have  given  theories  of  the  change  of  shape  due  to 

1  2 

the  effect  of  buoyancy,  notably  Penney-Price  and  Keller-Kolodner  .  These 
authors  represent  the  velocity  potential  cp  as  a  series  of  the  Legendre 
polynomials  Pq  (cos  9)  in  the  form: 

s  anPn^COS  9)R’^n+1^  (0.01) 

n=0  n  n 

where  R  is  the  distance  between  the  coordinate  origin  and  the  point  at 
which  the  velocity  potential  is  to  be  evaluated,  the  field  point.  Further, 
9  is  the  angle  between  the  vertical  axis  and  R,  and  the  &n’s  are  coeffi¬ 
cients  that  are  functions  of  time.  The  velocity  potential  can  be  inter¬ 
preted  as  the  sum  of  a  point  source,  a  dipole  and  higher  order  poles  all 
having  the  same  origin.  The  coordinate  origin  is  placed  on  the  axis  of 
symmetry  in  the  interior  of  the  bubble.  Since  the  potential  must  be 
evaluated  for  points  on  the  bubble  surface,  a  solution  is  possible  so 
long  as  the  bubble  surface  does  not  come  near  the  origin.  However,  the 
bubble  Jet  will  reach  every  point  in  the  interior  on  the  axis  of  symmetry 
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and  will  cause  Irremovable  singularities  as  R  goes  to  zero.  Therefore, 
the  representation  of  the  velocity  potential  falls  at  the  most  interesting 
and  important  part  of  the  motion  of  the  undervater  explosion  bubble. 

Initially,  the  bubble  expansion  is  radial  and  this  motion  is  known 
to  be  irrotational,  thus,  vortex  free.  However,  Snay  has  shown  that  as 
buoyancy  takes  effect  and  the  bubble  begins  to  migrate,  vorticity  i6 
generated. 

In  this  paper,  the  velocity  potential  will  be  represented  by  a  dis¬ 
tribution  of  point  sources  and  vortex  rings  on  the  surface  of  the  bubble, 
eliminating  singularities  from  the  interior  of  bubble  An  alternate 
approach,  as  suggested  by  Kelvin's  extent  ion  of  Green's  theorem,  is  a 
dj  stribution  of  point  sources  and  dipoles  over  the  bubble  surface,  but 
the  distribution  of  vortex  rings  is  simpler  in  the  analysi  i  if  the  vortex 
theorem  for  migrating  bubbles  is  applied. 

The  common  assumptions  of  bubble  theory  will  also  be  used  in  this 
study.  They  are: 

a)  The  bubble  motion  can  be  represented  by  a  velocity  potential  for 
points  outside  of  the  bubble  surface. 

b)  The  gas  in  the  bubble  expands  adiabatically  and  the  pressure  is 
uniform  throughout  the  bubble. 

c)  The  water  is  incompressible. 

In  this  paper,  the  water  is  assumed  to  be  unbounded. 

The  purpose  of  this  paper  is  not  to  solve  the  general  case  of  a 
migrating  bubble  which  changes  shape.  Rather,  the  effort  is  directed 
toward  deriving  the  equations  and  providing  the  mathematical  tools  for 
the  treatment  of  the  general  case  and  the  mathematical  treatment  of  the 
singularities  generated  by  the  distribution  of  point  sources  and  vortex 
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rings  on  the  surface  of  the  bubble,  ftie  resulting  equations  have  been 
tested  by  applying  them  tc  the  classical  case  of  a  migrating  spherical 
bubble. 
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CHAPTER  I 

BUBBLE  THEORY 

It  is  assumed  that  an  Incompressible,  unbounded  fluid  contains  a 
gas  bubble  within  it,  and  the  velocity  v^  of  the  fluid,  as  observed  from 
a  fixed  frame  of  reference  in  the  fluid,  is  assumed  to  be  derivable  from 
a  potential  function  cpf  which  satisfies  Laplace's  equation  so  that: 

Vf  =  -  vfcpf  (1.01) 

Vf\  =  o  (1.02) 

The  pressure  P  in  the  water  can  be  obtained  by  Bernoulli's  equation: 

Sj  - 1  (Vf >2  -  «ar  -  T0  -  f(t)  (1-03) 

where  the  subscript  f  references  the  spatial  and  time  coordinates  to  the 
fixed  reference  frame,  pQ  is  the  density  of  water,  g  is  the  acceleration 
due  to  gravity,  and  f(t)  is  the  integration  constant  arising  from  the 
integration  of  the  classical  hydrodynamic  equation  at  constant  time. 

The  positive  z  axis  is  assumed  to  point  vertically  upward.  It  is  further 
assumed  that: 


-•  0  as  -•  « 


(1.04) 


which  Implies  that: 


-»  0  and  VjflPf  -  0  as  Rf  -  »  (1.05) 
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where  is  the  distance  from  the  coordinate  origin  to  any  point  in  the 
fluid.  Statements  1.04  and  1.05  are  simply  the  mathematical  expressions 
for  the  assumption  that  the  fluid  is  undisturbed  at  Infinity  as  seen 
from  the  fixed  frame  of  reference,  consequently: 

P  -  PQ  as  R  -  «  (1.06) 

where  PQ  is  the  hydrostatic  pressure  of  the  fluid  for  2  =  0  at  infinity, 
thus: 


f(t) 


(1.07) 


The  Bernoulli  equation  can  now  be  written  as: 


Bcpf  1  ,  V2 

dtf  ’  2  ^  Vt^t '  “  gzf 


(1.06) 


The  pressure  P  is  continuous  across  the  boundary  of  the  bubble. 
Therefore,  the  pressure  in  the  bubble  must  be  the  same  as  that  in  the 
adjacent  water.  When  the  Bernoulli  equation  is  evaluated  at  the  bubble 
interface,  cui  interrelationship  between  the  pressure  P  in  the  bubble  and 
the  motion  of  the  surface  is  obtained. 

Since  the  bubble  stir  face  is  a  free  surface,  a  surface  that  moves  with 
the  fluid,  the  normal  velocity  of  the  bubble  surface  must  be  the  same  as 
that  of  the  adjacent  water.  According  to  Lamb^ ,  if  F  =  0  is  the  surface 
equation  of  the  free  surface,  the  surface  equation  must  satisfy  the  kine¬ 
matic  condition: 

VF’v  +  — ■  =  0  (l.09) 

ot 
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Die  unit  normal, 


of  any  surface  is: 


a  _  VF 
en  \VF\ 


(1.10) 


Therefore,  it  is  seen  that  if  expression  1.09  is  divided,  by  |vF|,  the 
normal  velocity  at  the  boundary  of  the  bubble  can  be  expressed  as: 


dF 

v  _  dt 
vn  '  ’  T?F[ 


(1.11) 


Assume  that  the  bubble  translates  in  the  positive  z  direction  vith 
a  velocity  B(t).  If  a  frame  of  reference  with  its  coordinate  origin  at 
the  bubble  center  is  rigidly  attached  to  the  bubble,  an  observer  at  the 
origin  would  see  a  stream  of  velocity  -  B  moving  past  the  bubble.  The 
relationships  between  the  variables  in  the  fixed  reference  frame  and 
those  in  the  moving  frame  are: 

zf  "  zo  +  B(t)  (1-12) 


tf  =  tm  (1.13) 

Vf  =  V,  (1.14) 


If  the  velocity  potential  in  the  fixed  reference  frame  is  expressed  in 
terms  of  the  variables  in  the  moving  reference  frame,  the  Bernoulli 
equation  can  be  written  as: 


m 


.  Bcp  1  _  2 

B  TT  '  I  (Vf  >  -  g(VB) 

m 


P-P 


(1.15) 
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where  the  superscript  m  shows  that  the  velocity  potential  in  the  fixed 
reference  system  is  expressed  in  terms  of  the  moving  coordinate  variables. 
The  relationship  between  the  velocity  in  the  fluid  as  seen  by  an  observer 
in  the  moving  frame  and  that  seen  by  an  observer  in  the  fixed  frame  is: 

-  vco  =  -  Fcp  -  Bk  (1.16) 

mm  nr  f 

From  expression  l.l6  it  is  easy  to  see  that: 

cp  =  cp“  +  Bz  (1*17) 

m  r  m 

where  the  subscript  m  refers  to  the  moving  reference  system.  Upon  making 
the  appropriate  substitutions,  the  Bernoulli  equation  expressed  in  terms 
of  the  velocity  potential  in  the  moving  frame  of  reference  and  the  moving 
frame  of  reference  and  the  moving  coordinates  is: 

H  -  BZ  +  |  B2  -  \  (Vcp f  -  g( z+B)  -  -^-2  =  o  (1.1B) 


The  Bernoulli  equation  can  be  written  as: 


1  *2 
+  iB 


1  2 
2  Vn 


1  2 

2  T 


P-P 

g(z+B)  -  —2  =  0 
po 


(1*19) 


where  v_  and  v  are  the  tangential  and  normal  velocities  as  seen  by  an 
T  n 

observer  in  the  moving  frame. 

In  order  to  solve  the  Bernoulli  equation,  an  expression  for  the 
velocity  potential  must  be  found. 

|i 

According  to  Lamb  ,  if  ro  and  cp'  are  any  two  single  valued  functions 
which  satisfy  Laplaces' s  equation  throughout  a  given  region,  then  using 
Kelvin's  extension  of  Green's  theorem: 


8 


(1.20) 
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where  5  Is  the  surface  of  a  boundary  enclosing  a  regi  n  r  the  fluid  and 
n  Is  the  coordinate  normal  to  S.  Let  v  be  the  velocity  potential  and 
cd'  =  l/R,  the  reciprocal  of  the  distance  of  any  point  of  the  fluid  to  a 
fixed  point  P  in  the  space  occupied  by  the  fluid.  It  is  necessary  to 
exclude  P  from  the  surface  of  integration  since  cp*  can  become  infinite  if 
P  is  included  in  the  region  of  integration.  This  may  be  done  by  describing 
a  small  sphere  about  P,  where  the  fixed  point  P  is  the  center.  If  £  refers 
to  the  surface  of  the  small  sphere  and  S  to  the  remaining  boundary,  expres¬ 
sion  1.20  gives: 


(1-21) 


2  2 

With  d(l/R)/dn  =  -  l/R  at  the  surface  £  and  d£  =  R  dv,  expression  1.21 
can  be  written  as: 


co  dv 


(1.22 ) 


As  R  -»  0  the  first  integral  on  the  right  hand  side  of  expression  1.22 
vanishes,  while  the  first  integral  on  the  left  hand  side  is  -  4ncpp,  where 
cop  is  the  value  of  cp  at  the  fixed  point  P.  Upon  rearranging  terms,  the 
velocity  potential  at  the  fixed  point  P  is: 


(1.23) 
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Expression  1.23  gives  the  value  of  cp  at  any  point  in  the  fluid  in  terms 
of  cp  and  dep/dn  at  the  boundary.  Since  the  velocity  v  is  defined  as: 

v  =  -Vcp  (1.24) 


the  normal  velocity  at  the  boundary  is: 


v 


n 


and  expression  1.23  can  be  written  as: 


.  l  jTvo 

"1 


®P  =  ^  * 


(1.25) 


(1.26) 


The  first  integral  represents  a  distribution  of  sources  and  the  second 
integral,  a  distribution  of  dipoles  on  the  boundary  S.  If  S  is  finite, 
both  integrals  vanish  as  R  -*  ®,  as  required  by  statements  1.04  and  1.05 
for  the  potential  in  the  fixed  system.  However,  cm  observer  in  the  moving 

•  A 

frame  of  reference  will  see  a  stream  with  velocity  -  Bk  flow  past  the  bub- 

•  A 

ble.  Therefore,  the  fluid  has  a  velocity  of  -  Bk  at  infinity  relative  to 
the  moving  frame  of  reference,  and  the  velocity  potential  of  the  stream  is: 

cp  ,  =  Bz  (1.27) 

^stream  '  ' 

The  velocity  potential,  of  the  bubble,  with  a  distribution  of  sources 
and  dipoles  on  the  bubble  surface,  as  seen  by  an  observer  in  the  moving 
frame  is: 


=  h 


(1.28) 


10 


NOLTR  66-211 


Maxwell"*  has  shown  that  any  magnetic  dipole  distribution  over  a  closed 
surface,  with  dipole  axes  directed  along  the  normals  of  the  surface,  may 
be  replaced  by  a  system  of  current-carrying  rings  distributed  over  the 
surface.  In  the  analogous  case  of  hydrodynamics,  the  distribution  of 
dipoles  can  be  replaced  by  a  system  of  vortex  rings  lying  in  the  surface. 
The  rings  have  a  common  axis  of  symnetry.  Such  a  representation  of  cp 
simplifies  the  analysis  when  the  vortex  theorem  for  migrating  bubbles  is 
applied. 

It  is  well  known  that  for  such  surface  distributions,  the  potential 
internal  to  the  boundary  of  the  surface  is  zero.  Therefore,  the  gradient 
of  the  potential  is  also  zero.  In  the  case  of  the  velocity  potential, 
the  velocity  internal  to  the  boundary  is  zero,  in  the  same  way  as  the 
electric  field  internal  to  a  charged  sphere. 

The  velocity  potential  for  a  single  vortex  ring  can  be  shown  to  be, 

6 

see  Milne-Thomson  : 


«v  =  $  (1,a9) 

where  r  is  the  circulation  and  Q  is  the  solid  angle  subtended  at  the 
field  point  by  the  diaphragm  which  is  enclosed  by  the  vortex  ring.  If 
y(s)  is  the  circulation  density  for  a  distribution  of  vortex  rings  on  the 
boundary  of  the  bubble,  then: 


cfy  =  ^  J  Y(s)nas  (1-30) 

where  the  integration  is  taken  over  the  external  boundary  of  the  bubble. 
The  circulation  is: 
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T  =  /  y(b )ds  (1.31) 

*0 

but: 

T  =jfvT(  s)ds  (1.32) 

where  v^,(s)  is  the  tangential  velocity. 

During  the  initial  expansion  of  the  bubble,  expressions  1.31  and 
1.32  are  zero  for  any  closed  curve  since  only  radial  flow  takes  place. 
However,  as  vorticity  is  being  generated  by  the  action  of  gravity,  the 
bubble  surface  becomes  a  surface  of  tangential  discontinuity  since  there 
is  a  translation  of  the  bubble  in  addition  to  the  radial  expansion.  If 
the  closed  path  is  totally  outside  the  bubble,  expression  1.32  is  zero 
at  all  times  since  the  fluid  motion  is  assumed  to  be  irrotational  at  all 
times.  Let  the  path  of  integration  be  on  the  external  and  internal 
boundary ,  thus : 


r  = 


(1.33) 


The  second  integral  must  be  zero  at  all  times  since  there  is  no  fluid 
motion  internal  to  the  bubble  boundary,  as  vas  previously  assumed.  UDon 
setting  expression  1.29  equal  to  1.31,  it  is  seen  that: 


( 1 . 3U ) 


since  the  path  of  integration  may  be  taken  as  small,  as  desired.  Therefore, 
the  potential  of  a  vortex  sheet  in  terms  of  the  tangential  velocity  at 
the  boundary  of  the  bubble  is: 
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<*V  =  ^  /vT(s)Wa  (1-35) 

•'s 

The  total  velocity  potential  due  to  the  source  sheet,  vortex  sheet  and 
otreaa  is: 


*? =  h 


v^sjftis  +  Bz 


(1.3b) 


The  time  derivative  of  expression  1.36  at  the  boundary  is  needed  for 
the  Bernoulli  equation. 

Let  the  velocity  potential  of  the  source  sheet  be  defined  as: 


(1-37) 


thus: 


*PS  _  1_  d_ 
dt  4tt  dt 


(1.38) 


Describe  a  small  sphere  about  the  point  P  on  the  boundary,  where  P  is  the 
center.  If  £  refers  to  the  surface  of  the  small  sphere  and  S  to  the  ori¬ 
ginal  boundary,  expression  I.38  gives: 


dt 


(1.39) 


If  d£  =  R  dv,  Expression  1.39  can  be  written  as: 


(1.1+0) 


13 


HOMR  66-211 


The  factor  dR/dt  in  the  second  integral  can  be  expressed  in  terms  of  the 
normal  velocity  as: 


dR 

dt 


v  sec  y 
n  1 


(l.4l) 


where  y  is  the  angle  between  the  normal  to  the  surface  S  at  the  point  P 
and  the  radius  vector  R.  Expression  1.40  can  now  be  written  as: 


.  i 
dt  5m 


Rdw  + 


1 

Im 


l 


v^sec  ydw 


1  iL 
5m  dt 


(1.42) 


The  term  sec  ydw  is  the  solid  angle  of  the  disk  subtended  at  a  point  on 
the  axis  of  symmetry  coinciding  with  the  normal  to  the  surface.  If  R  -•  0 
and  the  normal  is  the  outward  drawn  normal  to  the  surface,  the  first  inte- 

f 

gral  vanishes  and  J  sec  ydw  -♦  2n  and  one  obtains: 

I 


dt 


(1.43) 


The  time  derivative  of  the  velocity  potential  due  to  the  vortex  dis¬ 
tribution  will  be  treated  in  a  similar  manner  as  that  of  the  source  sheet. 
Describe  a  small  circle  in  the  surface  about  the  point  P,  where  P  is  the 
center  of  the  circle.  If  a  refers  to  this  circle  while  s  refers  to  the 
original  boundary,  the  time  derivative  of  the  velocity  potential  for  the 
vortex  distribution  can  be  written  as: 


dt  5Jn 


+  jL.  iL  f  vmftis 
4tt  dt  J  T* 

8-0 


(1.44) 


The  factor  dfl/dt  can  be  written  as: 
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BQ  _  BQ  BR 
Bt  "  BR  dt 


(1.45) 


and  dR/dt  can  be  expressed  in  terms  of  the  normal  velocity.  Expression 
1.44  can  nov  be  vritten  as: 

3T  =  In  030  +  55  ivTvnsec  Y  3R  40  *  55  It  1  vTnds  (lA6) 
a  a  s-a 

vhere  sec  y  BO/BR  is  the  projection  of  Bft/dR  onto  the  normal,  thus: 


an  _  bo 

8ec  Y  BR  B^ 


(1.47) 


vhere  n  is  the  coordinate  in  the  direction  of  the  normal.  Expression  1.46 
can  nov  be  vritten  as: 


Odrr  + 


f  dQ  ds 

I  V  V  —  - 

J  T  n  Bn  dn 

a 


dn  +  ~  rr 
4rr 


B_ 

Bt 


i 

;-n 


vTnds 


(1.46) 


vhere  ds/dn  is  the  inverse  of  the  slope  at  the  point  P.  If  do  =  rd0  and 
r  0,  the  first  integral  vanishes  and  the  time  derivative  of  the  velocity 
potential  for  the  vortex  distribution  is: 


B(X>y 

Bt 


1 

-  v  vm 
2  n  T 


cot  0  + 


1  —  f 

Bt  ✓ 


v  Qds 
T 


(1.49) 


vhere  the  angle  0  is  the  angle  betveen  the  normal  of  the  surface  and  the 
vertical  axis.  The  time  derivative  of  the  total  velocity  potential  at 
the  boundary  can  nov  be  vritten  as : 


dcpp 

Bt”" 


1  2 
2  Vn 


1 

2 


v  v„ 
n  T 


cot  0  + 


1 

Bt 


k 


Qds  +  Bz  (1.50) 
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Expression  1.50  refers  to  a  coordinate  system  moving  with  the  bubble 
and  having  its  origin  at  the  bubble  center.  Some  important  considerations 
concerning  the  coordinate  system  to  be  used  in  solving  the  problem  are: 

a)  The  problem  is  to  be  solved  in  terms  of  cylindrical  coordinates. 

b)  The  coordinates  of  any  point  on  the  surface  of  the  bubble  are  r, 
9,  and  z. 

c)  The  field  point  is  to  be  held  fixed  in  the  r-z  plane  and  its 

coordinates  are  r  and  z  - 

o  o 

d)  The  bubble  is  assumed  to  be  a  body  of  revolution  vith  its  axis 
of  symmetry  being  the  z  axis.  This  condition  implies  that  r 
and  z  are  independent  of  0. 
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CHAPTER  II 

VORTEX  THEORY  FOR  MIGRATING  BUBBLES 

The  Vortex-Theorem  for  migrating  bubbles  as  proposed  by  Snay  is  so 
vital  to  this  paper  that  it  will  be  derived  before  the  velocity  potential 
of  the  source  sheet  and  vortex  sheet. 

The  first  classical  hydrodynamic  equation  is  defined  as: 

♦  (v-v)v  =  a  -  ~  VP  (2.01) 

at  e  p 

where  v  Is  the  velocity  vector  in  the  fluid,  a  is  the  acceleration  due 

© 

to  any  external  forces  acting  on  the  bubble,  p  is  the  density  of  the 
fluid,  and  P  is  the  pressure  in  the  fluid.  To  an  observer  in  a  fixed 
coordinate  system  where  there  is  no  stream  present,  the  only  external 
force  acting  on  the  bubble  would  be  the  acceleration  due  to  gravity. 
However,  the  observer  in  a  moving  frame  observes  a  stream.  Therefore, 
the  acceleration  seen  by  an  observer  in  the  moving  frame  is  that  due  to 
the  stream  and  the  effect  of  gravity.  Expression  2.01  may  now  be  written 
as:  A 

+  (vxv)xv  =  (g+B)k  -  ~  v(v* v)  -  —  VP  (2.02) 

dt  2  p 

where  (v* v)v  has  been  expressed  in  terms  of  the  grad  v^  and  the  curl  v. 
If  expression  2.01  is  integrated  about  a  closed  path,  then: 

+  (vxv)xv  «ds  =^)(g+B)k*ds  -J)|^  v(v*v)  +  ^  vpj*ds  (2.03) 
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Since  the  fluid  is  incompressible,  the  last  integral  on  the  r.  h.  s.  can 
be  transformed  into  a  surface  integral  which  vanishes: 


JJj—  vxv(v*v)  +  -  VxVpJ*dS  =  0 


(2.04) 


Expression  2.03  becomes: 

+  (vxv)xvj*ds  =  j^(f5+B)k*ds  (2.05) 

If  the  flow  is  irrotational  everywhere,  as  in  the  case  of  radial 
flow  only,  expression  2.05  is  zero  no  matter  what  path  is  chosen.  However, 
if  there  exists  a  discontinuity  in  the  tangential  velocity  and  if  the  path 
of  integration  crosses  the  boundary  of  this  discontinuity,  expression  2.05 

is  different  from  zero.  As  was  discussed  in  Chapter  I,  the  velocity  of 

•  • 

the  fluid  internal  to  the  boundary  is  zero.  Therefore,  B  is  zero  internal 
to  the  boundary.  Uniformity  of  the  pressure  within  the  bubble  requires 
gravity  not  to  act  on  the  internal  fluid.  Thus: 


(g  ^B)k-ds  - 

c 

■'1 

where  the  x  .1  u.ion  is  over  the  external  boundary,  s  and  s.  are  the 
coord  i  r  1  wh-.M.  the  boundary  Is  crossed,  and  e  is  the  unit  tangential 

vector  at  the  boundary. 

A  A 

The  vector  (vxv)xv  is  zero  everywhere  except  for  the  infinitely  thin 
interior  of  the  vortex  sheet.  It  lies  in  the  plane  of  the  line  integral 
and  is  normal  to  the  boundary.  Hence: 


f 


A 


(g+B)k‘ 


ed: 

T 


(2.06) 


} 
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(2.07) 


where  e  is  the  unit  normal  to  the  vortex  sheet,  AT  is  the  thickness  of 
n 

the  sheet,  and  subscripts  1,  2  refer  to  the  points  of  crossing  on  the 
boundary.  Since  |0?xv)xv|  is  finite,  the  integrals  vanish  as  AT  -♦  0. 

The  integral  over  dv/dt  for  the  path  interior  to  the  bubble  is  zero 
since  the  velocity  is  everywhere  zero  internal  to  the  boundary.  The 
integral  of  dv/dt  over  the  paths  of  integration  interior  to  the  vortex 
sheet  can  be  written  as: 


lim 

AT-*0 


1 


lim 

AT-»0 


! 


at 


A  . 

e  dr 
n 


=  0 


2 


(2.08) 


where  the  subscripts  have  been  previously  defined.  Expression  2.05 
becomes : 


..  A 

J  (g+B)k- 


A 


ds 


(2.09) 


Vince  the  path  of  integration  can  be  made  as  small  as  desired,  the  inte¬ 
grands  in  expression  2.09  are  equal,  thus: 


St 


3. 


•  • 

+  B 


'/l+( r '  )d 


(2.10) 


where : 
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Bv  A 

at  ’  eT 


(ii.ll) 


A  A 

k*e„ 


(2.12) 


r'  =  (||)  (2-13) 

7 

To  illustrate  expression  2.10,  the  Herring-Zoller  equation  far  B  will 
be  derived  at  this  point.  The  time  rate  of  change  of  the  tangential 
velocity,  for  a  spherical  bubble,  at  the  boundary  is: 


^T  _  3  B  ,  3  (2.14 ) 

at  2  t - -  2  I - -r 

U+(r ')  a/l+Cr')^ 

where  a  is  the  normal  velocity  at  the  boundary  of  the  sphere  and  a  is  the 
radius  of  the  sphere.  Expression  2.10  for  a  spherical  bubble  can  be 
written  as: 

•  • 

•  •  TJn 

B  +  3  ~  =  2g  (2.15) 

3 

Upon  multiplying  both  sides  by  a  and  rearranging  terms,  expression  2.12 
can  be  put  into  the  form: 

(Ba3)’  =  2ga3  (2.16) 

Integrating  both  sides  of  expression  2.16  over  time  where  B  =  0  at  t  =  0 
gives: 


(2.17) 
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Expression  2.17  is  the  Herring-Zoller  equation. 

The  relationship  between  the  tangential  velocity  and  the  circulation 
density  of  the  vortex  rings,  established  in  Chapter  I,  leads  to  the  result 


ay(s) 

dt 


B  +  g 
>/l+(r  'f 


(2.15) 
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CHAPTER  III 

VELOCITY  POTENTIAL  OF  THE  SOURCE 
SHEET  AND  ITS  TIME  DERIVATIVE 


The  contribution  to  total  velocity  potential  by  the  source  sheet, 
expressed  in  cylindrical  coordinates,  can  be  written  as: 


(3.01) 


where  a  and  b  are  functions  o  time  and  are  the  value  of  z  at  the  top 
and  bottom  (poles)  of  the  bubble  where  r  =  0.  Expression  3«01  can  be 
written  in  the  form: 


(3-02) 


where  the  kernel  is  defined  as: 


X(r,z,ro,zo)  (3'°3) 

a  =  r"~  +  r  ^  +  (z  -z)^  (3-04) 

o  '  o 

0  =  ^rrQ  (3-05) 

According  to  Byrd^  ,  the  solution  to  an  elliptic  integral  of  the  type  in 


expression  3-03  Is: 
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(3.06) 


X(r,z,rQ,zo)  =  rpK(k) 

where  K(k)  is  the  complete  elliptic  integral  of  the  first  kind.  The 

modulus  k  and  the  variable  are  defined  as: 

k2  -  (3-07) 

u  =  {3’08) 

The  kernel  can  be  expressed  in  terms  of  k: 

K(r,z,ro,zQ)  =  K(k)k^  (3-09) 

o 

The  variable  r  is  defined  as : 

r  =  r(z,q1,q2,...,qn)  (3-10) 


where  the  q's  are  parameters  of  the  surface  and  are  functions  of  time 
only.  An  example  of  a  surface  parameter  would  be  the  radius  in  the  case 
of  a  sphere  or  the  semi-major  and  semi-minor  axis  in  the  case  of  an 
ellipse,  etc. 

The  partial  time  derivative  of  expression  3-02,  where  the  field 
points  are  fixed,  is: 


vn(r,z,t)K(r,z,ro, 


yi+o 


’)■ 


dz 


(3.11) 


+ 


K(r,z,rQ,zo 


r') 


2 


dz 
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& 

h  5f[*  Vr’z't)K(r’z’ro’zo^1+(r')‘:!|i  =  b 

-  a  vn(r,z)t)K(r,z,ro,zo)/l+(r');2  z  =  J 


where: 


ltK(r’z’Vzo) 


■  If  H[E(k)j1  +  s(-f)(£r)1 +  K(k)]kC 


(k')2  =  1  -  k2 


•  db 

b  =  dt 


a  = 


da 

dt 


(3.12) 


(3.13) 


(3.H) 


(3-15) 


and  k'  is  called  the  complementary  modulus.  If  k  and  k'  are  evaluated 
using  expressions  3.07  and  3-13,  the  positive  square  root  is  always  used. 
Expression  3«11  is  the  partial  time  derivative  of  the  velocity  potential 
for  the  source  sheet  at  any  point  in  the  fluid  excluding  points  on  the 
boundary  of  the  bubble.  As  was  shown  in  Chapter  I,  the  evaluation  of 
the  partial  time  derivative  of  the  velocity  potential,  at  the  boundary 
includes  a  l/2  v  ^(r  ,  z  ,t).  Therefore,  this  extra  term  must  be  added  to 
the  r.  h.  s.  of  express  ion  3«H  'when  the  partial  time  derivative  of  the 
velocity  potential  is  to  be  evaluated  at  the  boundary.  From  this  point 
on,  the  velocity  potential  of  the  source  sheet  and  its  time  derivative 
are  to  be  evaluated  at  the  boundary. 
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Recall  that  the  normal  velocity  is  derived  from  the  surface  equation 
F(r,z,t).  The  partial  operator  d/dt  operates  only  on  the  surface  para¬ 
meters.  Therefore,  the  normal  velocity  expressed  in  terms  of  the  time 
derivative  of  the  surface  parameters  is: 

vn  -  (3.16) 


vhere: 


|-F.(r,z,t) 

’  |VF(r,z,t)| 


(3. IT) 


The  partial  operator  d/dt  operating  on  r  and  VQ(r,z,t)  can  be  expressed 

as  ? 


dr  _  '  dr_ 

at  "  pi  dqi 


(3-18) 


d_ 

dt 


(3.19) 


It  must  be  remembered  that  d/dq^  operating  on  f  will  operate  on  r  keeping 
z  fixed,  but  the  partial  derivative  operates  on  F(r,z,t)  vhile  both  r  and 
z  are  fixed. 

The  velocity  potential  of  the  source  sheet  and  its  time  derivative 
can  now  be  expressed  as: 

<%(wt)  =  *Vo±  (3-20) 


d  t 

at 


r  ,z  ,t) 
o  o 


■  y + 


VW 


■♦a,  +0),  ) 


(3.21) 
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vher e : 


ao  =  Iff  /  fiK(z>z o)^/1+(r')2  dz  (3-22) 

i  */  a 


(3.23) 


(3.24) 

(3.25) 

(3.26) 

(3.27) 


The  new  notation  for  the  argument  of  the  kernels  shows  that  the  potential 
and  its  time  derivative  are  to  be  evaluated  on  the  boundary  where: 

ro  -  ro(VW  — >qn)  (3-S8) 


The  kernels  K(z,zq)  and  K^(z,zo)  have  singularities  at  z  =  zq.  As 
z  -»  z  ,  k  -*  1  and  K(k)  -*  «  while  s(k)  -*  1.  If  k  ~  1,  the  complete 
elliptic  integrals  K(k)  and  E(k)  can  be  represented  by  a  series  in  terms 
of  ln(4/k'),  Byrd'"^.  If  only  the  first  term  of  the  series  is  used  in 
each  case,  the  complete  elliptic  integrals  behave  as: 


K(k) 


(3.29) 
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E(k)  »  1 

Expand  r  in  terms  of  z  about  zq  such  that: 

r  =  r0  *  (r')0  (vz)  +  \  (r")02(vz) 

where  terms  higher  than  second  order  have  been  ignored  and: 


It  is  easily  shown  that  for  z  ~  z  ,  k'  behaves  as: 

o 


(3-30) 

(3.31) 

(3-32) 

(3-33) 

(3.31) 


Upon  substituting  the  r.  h.  s.  of  expression  3*3^  for  k'  in  expression 
3.29,  the  complete  elliptic  integral  K(k),  for  z  »  zq,  can  be  vriiten  as: 


(3.35) 


If  the  kernels  K(z,z  )  and  K^(z,zq)  are  expanded  in  terns  of  z  about 
and  terns  higher  than  the  zero  orde^  tern  are  ignored,  these  kernels 

behave  as: 
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V 


z,z 


o 


)  = 


(3.37) 


where  Y>(z,z^)  and  X^(z,zq)  are  the  approximate  kernels  that  represent  the 
behavior  of  the  original  kernels  about  the  singularity  and: 


(3-3B) 


(3.39) 


If  the  other  factors  in  the  integrals  of  expressions  3*22  -  3*25  are 
evaluated  at  z  =  zq  and  are  multiplied  by  the  appropriate  approximate 
kernel,  the  integration  of  the  approximate  integrand  over  z  is: 


(3.40) 

(3-41) 

(3-42) 

(3.43) 


where  a  ,  a.  ,  a-,  ,  and  Oi  are  approximate  coefficients  that  represent 
°i  2*  3i  4t 

the  behavior  of  the  original  coefficients  about  the  singularity  and: 


(3-44) 
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(3.^5) 


(3.46) 


If  the  approximate  coefficients  are  added  to  and  subtracted  from  the  ori- 
nal  coefficients,  the  singularity  is  removed  and  the  original  coefficients 
can  now  be  written  as: 


}  =^f  ["fiK(z,zo)/l+(r')2  -  f±  )i(z,zQ)fl+(r'  )0j<dz  +  aQ  (3-27) 

i  ^a  o  i 


2.  2TT 
i  'a 


J  K(z’z0^lt(r')2  '  (§^)0X(z'Z o)^I+(r')o‘J 


dz  +  o  (3-48) 


dz  +  o  (3.49) 
Ji 


cr, 


,  *  bj  W'-8. )/l+(r')2 '  fi  Ki(z’zo)4+(r’)02 

i  a  L  o  -1 


-  It/  fiK(*’zo>  £  fr1+(r')2)-  fi  *<2-zo>fc  /(-')a)} 

a  1  o  v  i  'o 


dz  (3.50) 


+  Ot, 


where  a  ,  a ,  ,  ,  and  a,  are  defined  as: 

°i  2i  3i  1 


t  -  -  ~  f-  /l+(r' )  d 
o,  2n  i  '  v  'o 
i  o 


!  (b-z  )/l+(r  * )/ 
(b.z  )ln  _ ^ - 2. 


(3.51) 


+  (z  -a)ln 
o 


(z  -a)/l+(r '  )  j2 


°  ^ ~ I  -  (b-a) 
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l+(r  * ) 


2 


(b-z  )in 
o 


(b-z  )/l+(r' )~ 


Ui 


(3-52) 


!(z  -a)/l+(r'  )J* 
+  (z  -a)ln - 


br 


-  (b-a) 


s  “  isr  fi  ^1+(rV(lj:) 

i  o  o  V  1i/o 


2(b-a)  + 


2(rr' ) 
x  'o 

l+(r')  ' - 


In 


z  -a 
o 


b-z 


(3.53) 


-  (b-zQ)ln 

(b-zo)/n-(r')oa 

(z  -a)ln 
0 

r 

iir 

0 

/  /  \ 

CVJ^ 

u 

+ 

r-l 

1  'r"* 

f'O  cr1 

1*? 

^||^J 

1 

II 

•H 

1  c7 
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The  coefficient  has  been  reserved  for  the  B  term  that  will  be  a 
1 

resalt  of  the  partial  time  derivative  of  the  potential  due  to  the  distri¬ 
bution  of  the  vortex  rings  on  the  surface  of  the  bubble.  The  reason  for 
this  type  of  notation  will  became  clearer  in  the  chapter  on  the  Bernoulli 


equation. 
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CHAPTER  IV 


VELOCITY  POTENTIAL  OF  THE  VORTEX  SHEET 
AND  TTS  TIME  DERIVATIVE 


The  contribution  to  the  total  velocity  potential  by  the  vortex  sheet, 
expressed  in  cylindrical  coordinates,  can  be  written  as: 


r , z, t )Q(r, z , rQ, 


dz 


(4.01) 


where  v?(r,z,t)  is  the  tangential  velocity  at  the  boundary  in  the  Moving 
frame  of  reference,  Q(r,z,ro,zQ)  is  the  solid  angle  of  the  vortex  ring 
subtended  at  the  field  point,  and  the  Halts  a  and  b  are  defined  in  the 
sane  manner  as  in  Chapter  III. 

The  partial  time  derivative  of  expression  4.01  is: 


z,t)Q(r,z,rQ, 


dz 


(4.02) 


As  was  shown  in  Chapter  I,  the  evaluation  of  the  partial  tiae  derivative 
of  the  velocity  potential  at  the  boundary  includes  an  additive  term  of 
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l/2  vn(ro,zQ,t)v^(ro,zo,t)  cot  0.  Therefore,  this  extra  term  must  be 
added  to  the  r.  h.  s.  of  expression  4.02  when  the  partied,  time  derivative 
of  the  velocity  potential  is  to  be  evaluated  at  the  boundary.  From  this 
point'  on,  the  velocity  potential  of  the  vortex  sheet  and  its  time  deriva¬ 
tive  are  to  be  evaluated  at  the  boundary. 

Recall  that  the  result  of  the  vortex  theorem  for  migrating  bubbles 
is: 

ft  yr.z.t)  -  llhr  <4-03) 

yi-Kr'r 


where  g  is  the  acceleration  due  to  gravity.  If  the  variable  r  is  defined 
in  the  same  manner  as  in  Chapter  III,  expression  4.02  can  be  written  as: 

ft  Vwt)  °  °i +  ^1^  Wvv0  cot  e)  (4-04) 


where  it  will  be  seen  later  in  the  chapter  that  Q(r,z,ro,zQ)  =  0  for  r  =  0 
and  this  makes  the  last  term  of  expression  4.02  zero  and: 


n 


1 


vT(r,z,t)n(z,zo)dz 


(4.05) 


a 


6. 

1 


|£.| -n(z.Zo)/T^7 


7.Z 


(4.06) 


dz 


(4. 07) 


The  solid  angle  q(z,zq)  is  the  same  as  in  expression  4.02,  but  the  new 
notation  for  its  arguement  shows  that  the  coefficients  are  to  be  evalu¬ 
ated  on  the  boundary  where: 
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ro  "  1'o(VVV"-’<!n)  (4-0B) 

The  solid  angle  of  a  vortex  ring  subtended  at  any  field  point  in  the 
fluid  is: 


n  ■ Js ,  fe(i?»)dS* 

where  R*  is  the  distance  from  the  field  point  to  any  point  in  the 
enclosed  by  the  vortex  ring  and  S*  is  the  area  of  the  plana  enclosed  by 
the  vortex  ring.  The  outward  drawn  normal  of  the  plane  is  the  unit  vector 

A 

k  in  the  positive  z  direction.  It  is  easily  shown  that: 


R*  =  ye"  +  r  2  +  (z  -z)“  -  2r  e  cos  0  (4.10) 


(vz) 


(z  -z)2  -  2r  e  cos 
v  o  o 


(4.11) 


where  e  is  the  variable  of  integration.  The  expression  for  the  solid 
angle  is  seen  to  be; 


n(r,z,ro,zo)  =  2  /  J 

0*/  o 


(z  -z)cilod0 


P 


+  r  ^  +  (z  -z)^  -  ;2er  cos  0 
O  0  o 


5 


w 


(4.12) 


vdiere  r  is  the  radius  of  the  vortex  ring  enclosing  the  plane.  Since  the 
integration  over  e  can  be  done  in  closed  form,  the  solid  angle  can  be 
expressed  as: 


33 


mum  66-211 


The  first  integral  can  be  integrated  over  0  in  closed  form  far  a  body  of 
revolution.  Therefore,  the  solid  angle  can  be  written  as: 


(vz> 

n(r,z,ro,zo)  -  ^ 

1  O 


-  a(vz)/ 


TT 


(i*.14) 


[ 


r  +  (z  -z )"'  -  r  r  c os  0  (10 
o  o  o  1 


'} 


r  2  .  v2  2  2.“|/~2 

r  +-  (z  -z )  -  r  cos  0  Yr 

_  o  o  o  _r  o 


2  \2 

+  r  +  (z  -z)  -  2rr  cos  0 

v  0  o 


After  using  the  method  of  partial  fractions  and  rearranging  terms,  the 
solid  angle  can  be  put  into  the  convient  fora: 


Q(r,z,ro,zo) 


=  ZTT 


±  1  - 


(z0'z) 


2n  •'  o 


J, 


X  +  r 


X  -  r 


X  +  r  cos  0 
o 


r  cos  9 
o 


-  9  cos  0 


-d0 


(4.15) 


where  a  and  B  are  the  same  as  in  expressions  3*04  and  3.O5  and  X  is 
defined  as: 


X  =  /ror"  +•  (zq-z)2  ( h .  16 ) 

The  integral  in  expression  4.16  can  be  broken  up  into  two  integrals  of  the 
type: 
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R(cos  6) 


S  cos  9 


d0 


where  R(coe  9)  is  a  rational  function  of  cos  6.  According  to  Byn?, 
the  integrals  can  be  reduced  to  Jacobian  elliptic  functions  of  the  form: 


0  “ 


.<2  2 
ak  sn  u 


(4.1b) 


■where  \j.  and  k  are  defined  in  expressions  3*07  and  3 .08.  If  the  substitu¬ 
tion  for  cos  0  in  expression  4.l8  is  made  in  the  tvo  integrals  resulting 
from  the  integral  in  4.15  and  the  indentities: 


nd2u 


1 

dn^u 


(4.19) 


2  ,  .2  2 

dn  u  =  1  -  k  sn  u 


(4.20) 


are  used,  the  solid  angle  can  be  written  as: 


fl(r,z,rQ,zo)  =  2n 


fi  1  -  (v2)kfc 

:)f 

o'"'  C 

K(k) 


dn^udu 

.  2  2 

1-0^  S’"  U 


(4.21) 


/X-r  \  dn^udu 

^~ro'^o  1-ou  ~sn“u 


where: 


(4.22) 


(4.23) 
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The  integrals  in  expression  4.22  are  complete  elliptic  integrals  of  the 

third  kind  and  have  five  possible  solutions  depending  on  the  value  of 

2  2 

and  •  The  condition  for  the  five  possible  solutions  are: 

I  o  <  -  ci^  <  00 

II  k2  <  <  1 

III  0  <  <  k2 

IV  1  <  a^2  <  ® 

V  complex  parameter 

2  2 

If  is  replaced  by  >  the  same  five  cases  exist  for  the  second 

2  2 

integral.  Case  V  is  excluded  since  and  cc,  811,6  always  positive  by 
definition.  It  is  easily  shown  that  always  satisfies  case  I  and 
will  satisfy  either  case  II,  III,  or  IV.  Therefore,  the  three  possible 
solutions  for  the  solid  angle  are: 


0  < 


2 

*1 


■-j 

<  k“ 


fl(r,z,r  ,z  )  =  2tt 


±  1  - 


(4.24) 


(a^-k^kjz^k). 

Yfi1^(l-o12)(k2-a12) 


B  = 


(4.25) 


(4.2 6) 
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k 


2 


<  1 


f}(r,z,ro,zo)  =  2n 


±  1  - 


(V2)k!/X+r 


2n,/r  r 

Y  o 


iX+r 


o/PPl 


K(k) 


(*.27) 


n(a1"-k"J)A0(^,k) 


2 )  J 


U-ro 


n(k“-a/)Ao(*,k) 


^a/(l-a./)(a/-k2) 


=  sin 


■1  “l 


2-k2 


a^k’)2 


(*.28) 


1  <  <  « 


fl(r,z,ro,zo)  =  L^rr 


±  1 


(Vz)k|A+r 


iX+r 


K(k) 


*L> 


(*.29) 


l(k‘--a12)K(k)z(n,k) 

/(aL"-l)(a1J-k2) 


A-r  / 

o/ 


n(k2-a22)A0((',k) 


°  -ir) 


If  the  partial  operator  d/dr  operates  on  the  solid  angle  in  expression 
4.14,  the  derivative  of  the  solid  angle  is  easily  shown  to  be: 


Q(r,z,r  , z  )  =  2(z  -z)  f 
dr  ’  o  o  o  y 


n 


rdQ 


c-0  cos  oj  J 


72 


(4.31) 


The  solution  to  the  integral  of  the  type  in  expression  4.32  has  been 
shown  by  Byrd^  to  be: 
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r 


a  9 


ja-3  cos  0^ 


375  “  A  E(k) 


C>.3<i) 


Upon  rearranging  terms,  expression  4-32  can  be  written  as: 


2_  n(r  -  r  ,  i  EfcL  -ii-fcal) 

ar  «r...r0.«0)  2^_  (kl)^r<)y 


(•’+•33) 


The  solid  angle  does  not  have  a  singularity,  but  the  derivative  of 

the  solid  angle  has  a  singularity  at  z  =  z  .  The  derivative  of  the  solid 

o 

angle  can  be  shown  to  behave  as: 


_ 2 _ 

(vz)[1+(r,)o"] 


(4.34) 


where  the  same  method  was  employed  here  as  in  Chapter  III  in  evaluating 
the  behavior  of  the  kernel  for  z  a  zq.  The  approximate  coefficient  for 

°61  18 : 


If  the  approximate  coefficient  is  added  to  and  subtracted  from  the  origi¬ 
nal  coefficient,  the  singularity  is  removed  and  the  coefficient  can  be^ 
written  as: 


a6*hJ  Vr’z>t)(i)f7  nU,zo)/i+(r')‘ 

i  a  \  1/ 


(4.36) 


v_(r 
T  o 


dz  +  a. 
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where: 


vm(r  ,zn,t) 


T 


-TT^/l+(r '  )  ‘ 


2  \*Vo 


In: 


z  -a 
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D-Z 


(7.37) 


39 


WO  I/PR  66-211 


CHAPTER  V 

TANGENTIAL  VELOCITY  OIL  THE  BOUNDARY  OF  A 
NONRIGID  BODY  OF  REVOLUTION  EXPANDING  AND 
MIGRATING  IN  A  LIQUID  MEDIUM 

The  tangential  velocity  at  the  boundary  of  a  migrating  and  oscil¬ 
lating  body  of  revolution  will  be  derived  in  the  moving  coordinate 
system. 

Assume  a  body  of  revolution  is  immersed  in  a  liquid  medium  and  the 
sin  fa ’O  S  of  the  body  divides  the  whole  space  into  two  regions  and  that 
both  regions  are  occupied  by  the  fluid.  If  the  fluid  external  to  the  body 
moves  past  the  body  with  a  velocity  B  in  the  positive  z  direction  and  the 
fluid  internal  to  the  body  is  at  rest,  the  surface  G  can  be  regarded  as 
a  surface  of 'discontinuity  of  the  velocity.  The  tangential  component  of 
the  velocity  can  be  produced  by  a  layer  of  vortex  rings  on  the  surface 
of  the  body. 

If  a  velocity  of  -  B  is  superpose!  onto  the  whole  flow,  so  that  the 
fluid  eternal  to  the  body  is  at  rest  at  infinity,  ther?  will  be  a  uniform 
parallel  flow  of  velocity  -  B  internal  to  the  body,  where  the  fluid  was  at 
rest  before.  This  flow  is  produced  by  the  vortex  layer  on  G.  The  problem 
is  to  determine  a  distribution  of  vortex  rings  on  the  surface  3  of  the 
body  which  produces  a  uniform  parallel  flow  internal  to  the  body  such 
that  the  tangential  velocity  at  the  field  point  P^ro,zQ)  on  the  boundary 
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It  has  been  shown,  by  Prandtl  ,  that  the  contribution  to  the  tangen¬ 
tial  velocity  by  the  vortex  on  which  the  field  point  is  located  is 
-  1/2  vT(ro,zo,t). 

The  induced  tangential  velocity  at  P(rQ,zo)  due  to  all  the  vortex 
rings,  except  the  one  located  at  the  field  point,  can  be  derived  through 
the  use  of  the  Biot-Survart  Lav: 


Hi 


A 

e. 


x  R 


v  (r,z,t) 


3 


J1  2 


1-4 

R 


1  e^ds^dSf 


(5-02) 


where : 

-  A  ,  ,  A  A 

R  =  i  (r  -r  cos  0  -  j  r  sin  0  +  k  (z  -z) 

o  iv  o 

A  A 

e^  =  -  i  sin  8  +  j  cos  0 
A  1  ^  )p  +  K 

/l+(r,)02 
ds1  =  / l+(r' )d  dz 


(5-03) 

(5.04) 


(5.05) 


(5.06) 


ds  =  rd0 


(5-07) 


Upon  substituting  expressions  5*03  -  5*07  into  expression  $.02  and  rear¬ 
ranging  terms,  the  induced  tangential  velocity  is  found  to  be: 


where  a  and  b  are  defined  in  the  same  way  as  in  Chapter  III  and: 
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X (r, z,r  ,  z  )  =  4  / 
T  o  o'  2J 


tt  r  -  r  cos  9 
o 


z  -z 

l--F-(r’) 


’^+r  +(z  -z)  -2rr  cos  9r' 

o  o  o  J 


rd9  (5.09) 


Expression  5-09  can  be  put  into  the  form: 


Vr,Z,ro,2o)  =  I 


7 


rr 


dQ 


o  (a-6  cos  e) 


37T-  v{1-¥!(r')oj/on(5.10) 


cos  9 


(a-9  cos  9) 


W 


d0 


where: 


2  x  2  ,  s2 
a  “  r  +  rQ  +  (zo-z) 


(5.11) 


and : 


8  =  2rrQ  (5.12) 

9 

According  to  Byrd  ,  the  solution  to  elliptic  integrals  of  the  type  in 
expression  5. 10  leads  to  the  solution: 


K_(r,z,r  ,z  )  =  77 
Tv  *  *  o’  o'  2 


z  -z 
o 


^  E(k)  +  r  r4l  -  -2—  (r’)  l  (5*13) 

La-3  o  (  r  of 


{■  ^  E(k)  ♦  *  K(k)  -  .4  E(k)} 


where: 


— —  (5-14) 

+  B 


Upon  expressing  a,  8,  and  yj  in  terms  of  r,  z,  r  ,  and  z  ,  expression  5*13 
can  he  expressed  as: 
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KT(r>z'ro'Zo)  = 


V(r  +r)2+(z -z)2  L 


K(k){l  -  ^r_!  (r-) 


•} 


(5-15) 


2r(r  -r)  z  -z  /  ^rr 

-  =<*>&• ♦  -  /.  7  -  -f-  (r'Ul  ♦  — rf—, 

(r  -r)  +(z  -z)  o  \  (r  -r)  +(z  -z) 


where  the  l/2  factor  is  now  taken  outside  of  the  integral  sign.  Expression 
5-15  can  he  simplified  further  to  give: 


yr,z,VZo)  =yf 


K(k){l 


Z  -Z  ) 

—  (r’)  J 

ro  °f 


(5.16) 


!,  /r  -r\/,  \2  z  -z 
,  .  If  o  Hk  i  o 

1  +  5l~jlk7  •  — 

\  o  /  '  '  o 


(r‘) 


1  + 

d\k' 


The  induced  tangential  velocity  is: 


vm  ( r  ,  z  ,  t )  - 
T  '  o  o’  ' 


^r/l+Cr1)/ 


/  vT(r,z,t)J^(r,z,r  ,z  )/l+(r‘  )'*  dz  (5. IT) 


where  the  kernel  J^,(r,z,ro,zo)  is  that  of  expression  5-1 6. 

If  the  tangential  velocity  at  P(ro,zQ)  due  to  the  fluid  flow  internal 
to  the.  boundary  is  set  equal  to  the  sum  of  -  l/2  v^(rQ,zo,t)  and  the  induced 
tangential  velocity,  the  tangential  velocity  v^r^z^t)  is  found  to  be: 


(5.1«) 
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The  kernel  in  the  integral  of  expression  5.13  has  a  singular  point 
at  z  =  zq.  If  the  same  method  is  used  here,  as  vas  employed  in  Chapter 
III,  to  remove  the  singularity,  the  approximate  kernel  is  found  to  be: 


uhere: 


(5*20) 


The  tangential  velocity  at  the  boundary  with  the  singularity  removed  is: 
vT(ro,zQ,t)  =  2B  -“J*  | vT( r , z , t  )K( r ,  z , rQ , zq  )/l+( r '  )2  ( 5 - 21 ) 


-  vT(VVt)KT(r'a»Vzo)Vl+(r,)o2  dz 


z  -  ■  ■  -1 - + 

.  )ll+(r’ )  2 


VT(ro,Zot) 


,  „  flVaH1+<rV).  (l^0IVlHr')o 

(zo-a)lnj - 5^ - J+  (b-zo)ln| - - 


»  \  2 


lb“z  lY1+(r') 


•  \  2 
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CHAPTER  VI 

THE  BERNOULLI  'EQUATION 


Recall  that  the  Bernoulli  equation  in  the  moving  coordinate  system 
was  written  as: 


It  ctl(ro’zo’t)  •  \  *  \  •  I  vn2(ro’*o’t)  +  vT2(ro'zo’t)j  (6- 

L 


.01) 


1  p-p 

2  +  B(t) - 2  =  0 


If  the  time  derivative  of  the  velocity  potentials  for  the  source  sheet 
and  the  vortex  sheet  are  added  together  to  give  the  time  rate  of  change 
of  the  total  velocity  potential  and  vn(ro,zQt)  is  expressed  as  a  function 
of  the  q  ,  the  Bernoulli  equation  can  be  written  as: 


q.o  +  q ,2(o0  +  a?  +  Ok  )  +  q,(0c  +  Oc  +  07  + 

1  °i  1  *i  Hi  ?i  i  'i 


(6.02) 


I  fl  ¥T(Wt)cot  9o)] 

O  -J 


p-p 


|  B2  -  |  \(r0,z0,t)  *  0l  -  g[z0  ♦  B(t)l  -  *  o 

J  o 


The  Bernoulli  equation  must  be  evaluated  for  i  field  points  on  the  surface 
of  the  bubble  since  a  solution  must  be  found  for  each  q^.  The  system  of 
equations  can  be  put  into  the  form: 
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where: 


1*2  * 
q,  (a0  +  a,  +  ok  )  '■<L{  (O'-  ♦  Os  o7 

.  *)i  ^  5ji  ji  ji 


(6.0>  ) 


2  fJi0VT(roJ’ZoJ’t)cot  9o^]  +  5  ^ 


-  I  ’,T2(roJ’ZoJ't)  +  alj  ‘  8[v,  +  B(t)]  *  (~pf)  =  ° 


and  the  subscript  j  Is  for  the  evaluation  of  the  equation  at  the  j'th 
field  point.  Let  X  be  the  variable  matrix,  C  the  coefficient  matrix  and 
Q  the  matrix  whose  elements  are  the  Q^,  the  system  of  equations  can  now 
be  written  as: 


CX  +  Q  =  0 


(6.05) 


where: 


(6.06) 

(6.07) 


(6.08) 


(6.09) 
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If  the  inverse  of  the  coefficient  operates  on  expression  6.08,  the  equation 
is  transformed  into: 

X  =  -  C'^Q  (6.10) 

Since: 

C-1C  =  I  (6.11) 

and: 

IX  =  X  (6.12) 

vhere  I  is  the  unitary  matrix.  The  resulting  matrix  on  the  r.  h.  s.  of 
expression  6.14  is  the  solution  matrix  v&ere  each  element  in  the  matrix 
is  the  solution  to  the  corresponding  element  in  the  X  matrix.  Expression 
6.14  can  now  be  written  as: 

X  =  S  (6.13) 

vhere: 

S  =  -  C~\  (6.14) 

(6.15) 

and: 

(6.16) 
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Therefore,  the  solution  for  the  q's  can  be  arrived  at  in  the  same  manner 
as  solving  a  set  of  linear  algebraic  equations.  Once  the  solution  of  the 
q's  is  known,  the  q^'s  and  q^'s  can  be  solved  for  by  integrating  the 
variables  over  time. 

Now  that  the  qi's  and  q^s  are  known,  the  shape  of  the  bubble,  the 
normal  velocity,  and  the  tangential  velocity  at  the  boundary  can  be  deter¬ 
mined.  Therefore,  the  velocity  potential,  the  time  derivative  of  the 
velocity  potential,  and  the  spatial  derivatives  of  the  velocity  potential 
can  be  found  at  any  point  in  the  fluid.  The  Bernoulli  equation  can  now 
be  solved  for  the  pressure  at  any  point  in  the  fluid. 
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CHAPTER  VII 

THE  MIGRATING  SPHERICAL  BUBBLE 


The  surface  equation  for  a  sphere  in  the  moving  reference  frame  is: 

F(r,z,t)  =  r2  +  z2  -  a2  =  0  (7.01) 

vhere  a  =  a(t)  is  the  radius  of  the  sphere.  Since  there  is  only  one 
surface  parameter  =  a  (compare  with  expression  3«10)>the  normal  velo¬ 
city  at  the  boundary  of  the  sphere  is: 

vQ  =  a  (7.02) 

and  the  function  defined  by  expression  3*17  becomes: 

f 1  =  1  (7.03) 


The  velocity  potential  for  a  migrating  spherical  bubble  in  cylindrical 
12 

coordinates  is  : 


cp(wt}  =  A 


)/r  2+z  2 

»  n  c\ 


a3z 


(r  c+z  ) 
'  o  o 


^172 


+  Bz 


(7.04) 


where  r  and  z  are  the  field  point  coordinates  and  B  =  B(t)  is  the  velo- 
°  o 

city  of  migration.  The  tangential  velocity  at  the  boundary  of  a  moving 
sphere  is: 


(7.05) 
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The  time  derivative  of  the  velocity  potential  due  to  the  source  and 
vortex  sheets  is  obtained  from  expressions  3*21  and  4.04: 

aacQ  +  Bac^  +  a^c2  +  Bac^  +  Bzq  (7*06) 

where:* 


c 

o 


(7.07) 


C1  =  °1  +  °7 


(7.08) 


c2  =  a3  +  %  +  a5  +  2 


(7.09) 


c 


3 


(7.10) 


The  are  given  by  the  expressions  3*22  through  3.25  and  4.05  through 

4.07.  The  l/2  term  in  c^  comes  from  the  additive  term  included  in  the 
time  derivative  of  the  velocity  potential  of  the  source  sheet  and  the 
3A(zQ/a)  term  is  from  the  additive  term  included  in  the  time  derivative 
of  the  velocity  potential  of  the  vortex  sheet.  The  normal  and  tangential 
velocities  have  been  substituted  into  the  ai  coefficients  and  the  a,  B, 
and  B  variables  have  been  factored  out  to  give  cq,  c^,  c2,  and  c^.  If 
the  integrands  of  <3^  and  are  combined  as  veil  as  those  of  and 
the  coefficients  in  expression  7.07  can  be  written  as: 


c 

o 


1 

2na 


.  z 

j  °<v°o)dz + n°o-®o)dz + it 

-n  ^  7. 


(7.11) 


*  The  ci  coefficients  correspond  to  the  aQ  and  an  coefficients  resulting 
from  the  partial  time  derivative  of  expression  0.01  as  follovs:  c  =  a  , 
ci  a  ax,  c2  =  aQ,  and  c  - 
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C1  =  5na 


i 


o  fa 

Q(zQ,z)dz  +  J  fi(z  ,z)dz 
-a  •'z 

o 


(7.12) 


c2  2n 


*z 

t  o  f 

J  (G2-G2)dz  +  j  (G2-G2)dz  +  Ig 
.  “a  Zo 


(7.13) 


+ 


1 

2 


_  3 

c3  '  ^ 


z 

j  °(G3-G3)d2  +  J  (G3-G3)dz  +  I. 


(7.H.) 


where  Gq,  G2>  and  G^  are  the  original  integrands  with  singularities  at 

z  =  zq;  Gq,  G  ,  and  G^  are  the  approximate  integrands  far  z  a  zqJ  Iq,  I2, 

and  I_  are  the  analytic  integrals  and  fi(z  ,z)  is  the  solid  angle  of  a 
3  o 

vortex  ring  at  z  subtended  at  the  field  point  z  .  The  difference 

o 

between  the  original  integrand  and  the  approximate  integrand  is  called 
the  compensated  integrand.  The  original  integrands,  approximate  integrands, 
and  the  analytic  integrals  for  a  spherical  migrating  bubble  are: 
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(7-22) 


(7.23) 


If  expression  7.0U  is  differentiated  with  respect  to  time,  keeping 
the  field  point  variables  fixed,  and  evaluated  on  the  boundary,  the 
coefficients  cq,  c^,  c^,  and  become: 

=0  =  1  (7.21*) 
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,  ■  i  0 


(7.25) 


c2  =  2 


(7.26) 


■m 


(7.27) 


The  tangential  velocity  at  the  boundary  of  a  spherical  migrating 
bubble  is: 


VT  =  BcT 


(7.28) 


vhere  the  tangential  velocity  in  expression  7-05  has  been  substituted 
into  expression  5.21  and  B  has  been  factored  out.  The  coefficient  c^,  is 
written  as : 


where : 


(7.29) 


(7.30) 
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lhe  tangential  velocity  obtained  from  the  velocity  potential  in  expression 
7.04  evaluated  at  the  boundary  shows  the  coefficient  to  be: 


”  2 


(7.33) 


The  coefficients  in  expressions  7.11  through  7*15  and  expression 
7.29  were  programmed  by  this  writter  for  the  IBM  7090  computer  and  evalu¬ 
ated  for  several  field  points  on  the  boundary  of  a  unit  sphere.  Since 

all  (G  -  G)  integrands  are  continuous  functions  of  z  such  that  they  are 
zero  or  finite  at  z  -  zq,  the  integral  can  be  obtained  by  a  numerical 

quadrature;  e.g. ,  Gauss'  formula.  Hie  Gaussian  formula  has  an  advantage 

in  that  the  integrand  is  not  evaluated  at  the  limits  of  integration.  Using 

zq  for  one  of  the  limits,  the  integrand  is  computed  for  values  of  z  /  zq, 

thus  eliminating  the  need  of  determining  the  integrand  for  z  =  zq.  If 

the  integrand  is  evaluated  at  n  points,  Gauss'  formula  integrates  a 

(2n  -  l)'th  polynomial  exactly. 

Table  la  is  a  list  of  computed  values  for  cq  and  cg,  vhere  an  8-point 
Gaussian  quadrature  was  used  for  each  integration  range.  The  accuracy  of 
the  computed  coefficients  cq  and  c^  decreases  somewhat  as  the  field  point 
moves  away  from  the  equator  of  the  sphere  where  zq  =  0,  The  reason  for  this 
loss  of  accuracy  can  be  qualitatively  seen  in  figures  3  and  4.  These  are 
plots  of  the  original  integrand  Gq,  approximate  integrand  Gq,  and  compen¬ 
sated  integrand  (Gq  -  Gq)  as  a  function  of  the  integration  variable  z. 

The  curve  of  the  compensated  integrand  in  figure  4  (zq  =  0)  appears  to 
have  a  parabolic  shape  which  the  Gaussian  formula  would  integrate  exactly. 
However,  the  curve  of  the  compensated  integrand  in  figure  3  (zq  =  0*5) 
appears  to  have  higher  inclinations  in  the  region  of  z^  than  are  commen¬ 
surate  with  a  polynomial  and  a  lower  accuracy  of  integration  is  to  be 
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Table  lb  is  a  list  of  the  computed  values  of  and  as  compared 

with  the  true  values  calculated  using  expressions  7*25  and  7 .27.  Figures 

5  and  6  are  plots  of  the  integrands  Gy  Gy  and  (G^  -  G^).  Judging  from 

the  shape  of  the  curves  for  r.tae  compensated  integrands,  the  accuracy  of 

c_  would  be  expected  to  be  Is  so  for  z  =0.5  than  for  z  =0.  Hie  numeri- 
3  oo 

cal  results  for  Cy  c ,  and  show  the  same  trend  of  accuracy. 

Although  the  accuracy  of  the  coefficients  is  sufficient  for  practical 
purposes,  it  could  be  improved  as  follows :  Figures  3-6  suggest  employment 
of  four  regions  of  integration,  thus: 


c 


i 


-e 

G^dz  + 


,z 

j  (G1-Gi)dz 


(7.30 


+ 


)dz  + 


J\  +  *o. 

z  +e  i 

o 


where  an  e  is  chosen  such  that  the  region  of  high  inclination  is  covered 
by  the  two  inner  integrals.  The  compensated  integrand  is  not  necessary 
for  the  two  outer  regions  of  integration  since  there  are  no  singular  points. 
This  method  has  not  been  tried  yet. 

The  c2  coefficient,  like  the  coefficient,  includes  an  additive 
term  when  evaluated  on  the  boundary  of  the  bubble  (1.43).  However,  if 
the  coefficient  is  evaluated  at  points  off  the  boundary,  the  additive 
term  is  not  needed.  Also,  the  method  of  compensating  integrands  is  not 
needed  since  there  are  no  singularities  for  field  points  off  the  boundary. 
Figure  7  is  a  plot  of  the  integrand  G^  as  a  function  of  the  integration 
variable  z,  for  zq  =  0.  Hie  distance  of  the  field  point  from  the  boundary 
is  rQ  -  1.  Since  G^  is  symmetrical  about  the  G2  axis,  only  values  for 
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z  i  0  are  shown.  The  value  of  G„  is  finite  for  z  =  0  and  r  >  1.  The 

d  Q 

integral  of  G^  yields  the  value  of  c 2  shown  in  figure  8.  The  value  of 

is  extrapolated  to  2  at  rQ  =  1  which  is  the  correct  value  of  c  at  the 

boundary.  As  rQ  -•  1,  the  curve  becomes  more  and  more  pointed  and  the 

integral  of  G„  increases  so  that:  c„  -•  2  for  r  =  1  +  e.  However, 

d  a  o 

for  rQ  =  1,  the  G  curve  lies  very  close  to  the  axes  in  figure  7  and  the 
integral  of  G 2  decreases  in  value  such  that  cg  =  3/2.  The  additive  term 
of  l/2  gives  the  correct  value  of  2. 

CONCLUSIONS 

The  above  numerical  analysis  pertaining  to  a  spherical  bubble  demon¬ 
strates  that  an  underwater  explosion  bubble  can  be  represented  by  a  dis¬ 
tribution  of  sources  and  vortex  rings  on  the  surface.  Since  the  equations 
were  derived  v'thout  regard  to  the  shape  of  the  bubble,  they  should  be 
applicable  to  the  case  of  a  nonspherlcal  bubble.  This  method  would  not 
encounter  the  problems  of  the  model  consisting  of  a  single  point  source 
and  a  single  dipole  located  on  the  axis  of  symmetry  which  interferes  with 
the  Jet  rising  into  the  interior  of  the  bubble. 

It  is  planned  in  the  near  future  to  apply  the  theory  and  mathematical 
tools  developed  in  this  paper  to  the  case  of  a  nonspherlcal  bubble  with 
a  Jet.  The  differential  equations  for  the  q's  are  to  be  solved  using  the 
Adams -Moulton  predictor-corrector  method  or  that  of  Runge-Kutta  .  A 
generalized  program  using  these  two  methods  has  been  developed  at  the 
Naval  Ordnance  Laboratory  and  has  been  used  successfully  in  the  past. 

The  model  of  a  Jet-forming  bubble  would  allow  for  a  detailed  and 
complete  study  of  bubble  phenomena  associated  with  underwater  explosions. 

A  study  of  this  type  would  greatly  enhance  the  present  knowledge. 
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TABLE  la 

COMPUTED  VALUES  OP  c  AND  c  OBTAINED  FROM  EXPRESSIONS  7. 11 

o  2 

AND  7.13  EVALUATED  ON  THE  BOUNDARY  OF  A  UNIT  SPHERE 

These  coefficients  are  needed  for  the  time  derivative  of  the  velocity 
potential  of  a  migrating  spherical  bubble  and  are  associated  with  the 
potential  due  to  a  source  sheet.  The  true  values  of  cq  and  c  are  1  and 
2  as  obtained  from  expressions  7*24  and  7.2 6. 


Field  Point 

Coordinate  z 

0 

c 

0 

Computed  Values  Of: 

C2 

0.9 

0.99980 

2.00251 

0.8 

0.99995 

2.00058 

0.7 

0.99998 

2.00036 

0.6 

0.99999 

2.00010 

0.5 

0.99999 

2.00001 

0.4 

0.99999 

1.99998 

0.3 

0.99999 

2.00000 

0.2 

0.99999 

2.00000 

0.1 

0.99999 

2.00000 

0.0 

1.00000 

2.00000 

-0.1 

0.99999 

2.00000 

-0.2 

0.99999 

2.00000 

-0.3 

0.99999 

2.00000 

-0.4 

0.99999 

1.99993 

-0.5 

0.99999 

2.00001 

-0.6 

0.99999 

2.00010 

-0.7 

0.99998 

2.00036 

-0.8 

0.99995 

2.00058 

-0.9 

0.99980 

2.00251 
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TABLE  lb 

COMPUTED  VALUES  OF  c±  AND  c  OBTAINED  FROM  EXPRESSIONS  7.12 
AND  7.14  EVALUATED  ON  THE  BOUNDARY  OF  A  UNIT  SPHERE 

These  coefficients  are  needed  for  the  time  derivative  of  the  velocity 
potential  of  a  migrating  spherical  bubble  and  are  associated  vith  the 
velocity  potential  due  to  the  vortex  sheet.  Tlie  true  values  of  c^  and 
are  obtained  from  expressions  7.25  and  7.27* 


Field  Point 
Coordinate  z 

o 

Computed  Value 
of  c^ 

True  Value 
of 

Computed  Value 
°f  c3 

True  Value 
of 

0.9 

0.44648 

0.45 

1.35006 

1.35 

0.8 

0.39382 

0.40 

1.19985 

1.20 

0.7 

0.34142 

0.35 

1.04972 

1.05 

0.6 

0.29440 

0.30 

0.89989 

0.90 

0.5 

0.24353 

0.25 

0.75000 

0.75 

0.4 

0.19555 

0.20 

0.60009 

0.60 

0.3 

0.14701 

0.15 

0.44997 

0.45 

0.2 

0.09941 

0.10 

0.29995 

0.30 

0.1 

0.05008 

0.05 

0.14998 

0.15 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.1 

-O.O5OC8 

-0.05 

-0.14998 

-0.15 

-0.2 

-0.09941 

-0.10 

-0.29995 

-0.30 

-0.3 

-0.14701 

-0.15 

-0.44997 

-0.45 

-0.4 

-0.19555 

-0.20 

-0.60009 

-0.60 

-0.5 

-0.24353 

-0.25 

-0.75000 

-0.75 

-0.6 

-0.29440 

-0.30 

-0.89989 

-0.90 

-0.7 

-0.34142 

-0.35 

-1.04972 

-1.05 

-0.8 

-0.39382 

-0.40 

-1.19985 

-1.20 

-0.9 

-0.44648 

-0.45 

-1.35006 

-1.35 
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TABLE  2 


COMPUTED  VALUES  OF  c^,  OBTAINED  FROM  EXPRESSION  7.29 
EVALUATED  ON  THE  BOUNDARY  OF  A  UNIT  SPHERE 


This  coefficient  is  needed  for  the  tangential  velocity  of  a  migrating 
spherical  bubble.  The  true  values  of  c^,  are  obtained  from  expression  7.33* 


Field  Point 
Coordinate  z 

0 

Computed  Value 

cT 

True  Value 

of 

0.9 

0.65373 

0.65383 

0.8 

0.89996 

0.90000 

0.7 

1.07119 

1.07121 

0.6 

1.19999 

1.20000 

0.5 

1.29903 

1.29904 

0.4 

1.37477 

1.37477 

0.3 

1.43091 

1.43091 

0.2 

1.46969 

1.46969 

0.1 

1.49248 

1.49248 

0.0 

1.50000 

1.50000 

-0.1 

1.49248 

1.49248 

-0.2 

1,46969 

1.46969 

-0.3 

1.43091 

1.43091 

-0.4 

1.37477 

1.37477 

-0.5 

1.29903 

1.29904 

-0.6 

1.19999 

1.20000 

-0.7 

1.07119 

1.07121 

-0.8 

0.89996 

0.90000 

-0.9 

0.65373 

0.65383 

59 


NOLTR  66-211 
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Figure  3  THE  INTEGRANDS  Gq,  Gq,  AND  (Gq  -  Gq)  IN  EXPRESSION 
7.11  PLOTTED  AS  A  FUNCTION  OF  THE  INTEGRATION  VARIABLE  z 

'The  integrands  are  evaluated  at  the  boundary  of  unit  sphere  for 
zo  =  °*5* 
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Figure  4 

Same  as  figure  3  for  z  =  0. 

o 
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Figure  5  THE  INTEGRANDS  G^,  G^,  AND  (G^  -  G^)  IN  EXPRESSION 
7.14  PLOTTED  AS  A  FUNCTION  OV  THE  INTEGRATION  VARIABLE  z 

The  integrands  are  evaluated  at  the  boundary  of  a  unit  sphere  for 

0. 


) 
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Figure  6 

Game  as  figure  5  for  zq  -  0.5* 
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